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On the degree conjecture for separability of 
multipartite quantum states 



Ali Saif M. Hassan and Pramod S. Joa^ 
Department of Physics, University of Pune, Pune, India-411007. 

We settle the so-called degree conjecture for the separability of multipartite quantum 
states, which are normalized graph Laplacians, first given by Braunstein et al. [Phys. 
O ' Rev. A 73, 012320 (2006)]. The conjecture states that a multipartite quantum state is 

^ ■ separable if and only if the degree matrix of the graph associated with the state is equal 

. to the degree matrix of the partial transpose of this graph. We call this statement to be 

[T I \ the strong form of the conjecture. In its weak version, the conjecture requires only the 

^ \ necessity, that is, if the state is separable, the corresponding degree matrices match. We 

prove the strong form of the conjecture for pure multipartite quantum states, using the 
modified tensor product of graphs defined in [J. Phys. A: Math. Theor. 40, 10251 (2007)], 



O ■ as both necessary and sufficient condition for separability. Based on this proof, we give 

^ ! a polynomial-time algorithm for completely factorizing any pure multipartite quantum 
state. By polynomial-time algorithm we mean that the execution time of this algorithm 

^ ■ increases as a polynomial in m, where m is the number of parts of the quantum system. 

Q ■ We give a counter-example to show that the conjecture fails, in general, even in its weak 

^ . form, for multipartite mixed states. Finally, we prove this conjecture, in its weak form, 

Oh! for a class of multipartite mixed states, giving only a necessary condition for separability. 
> 

c3 : PACS numbers:03.67.-a,03.67.Mn 



I. INTRODUCTION 



The problem of detection and quantification of entanglement of multipartite quantum 
\ states is fundamental to the whole field of quantum information and in general to the 

physics of multicomponent quantum systems. Whereas entanglement of pure bipartite 
states is well understood, the classification of mixed states according to degree and char- 
acter of their entanglement is still a matter of intense research [3,4]. Presently, the most 
successful approach for bipartite and multipartite entanglement is via PPT criterion [5], 
CCNR criterion [6], positive maps and entanglement witnesses [7,8], Bloch representation 
of quantum states [9,10], covariance matrices [11] and local uncertainty relations [12]. All 
these approaches lead to either necessary or sufficient (but not both) criteria for separa- 
bility. Recently, a combinatorial approach to the separability problem is initiated by S. 
L. Braunstein, S. Ghosh and S. Severini [13] and extended by us [2]. In this approach the 
density matrices are coded in terms of graphs with the idea of using the graph topology 
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and operations on graphs for detection, classification and quantification of entanglement 
in the multipartite quantum states. S. L. Braunstein, S. Ghosh, T. Mansour, S. Severini, 

and R. C. Wilson [1] have made a conjecture, called degree conjecture, for the separability 
of multipartite quantum states. The conjecture states that a multipartite quantum state 
is separable if and only if the degree matrix of the graph associated with the state is equal 
to the degree matrix of the partial transpose (with respect to a subsystem, see below) 
of this graph. We call this to be the strong form of the conjecture. In its weak version 
it requires only the necessity, that is, if the state is separable the corresponding degree 
matrices match. We prove the strong version for a m-partite pure state (section 11) and 
give a polynomial-time algorithm, for factorization of a m-partite pure state (section 111). 
We show that the conjecture fails, in general, for mixed states (section IV). Finally, we 
prove the weak version of the conjecture for a class of multipartite mixed states (section 
V). 



II. THE SEPARABILITY CRITERION AND ITS PROOF 



For the definition of a weighted graph with real or complex weights, denoted {G,a), 
we refer the reader to ref.[2]. We denote by V{G,a) and E{G,a) the vertex set and the 
edge set respectively of a weighted graph. The degree of a vertex v e V{G, a) is denoted 
by D„ and 5{G,a) = is the degree sum of the graph. The degree matrix of a weighted 

graph is denoted by A{G,a) [2]. The combinatorial Laplacian of the weighted graph is 
denoted L{G,a) and the generalized Laplacian of (G, a) is denoted Q{G,a) [2]. If the 
generalized Laplacian of a graph (G, a) is positive semidefinite we can define the density 
matrix of the graph (G, a) as o"(G, a) = — Q{G, a). Conversely, given a density matrix 

we can assign a graph to it [2]. A real weighted graph gets assigned to a density matrix 
with all real elements, while a complex weighted graph is assigned to a density matrix 
with complex off-diagonal elements. The vertices of the graph are labeled by the elements 
of the standard basis in the state space of the multipartite quantum system which is used 
to set up the density matrix. Also, in what follows we use the definition and properties 
of the modified tensor product of graphs proved in ref . [2] . 

Let (G, a) be the graph corresponding to a m-partite pure state in d — did2 ■ • ■ dm 
dimensional Hilbert space H — T-f^^ ® H^'^ • • • (8) Ti'^'" , where di is the dimension of 
Tl'^\ 1 < i < m. Each vertex of (G, a) is labeled by an m tuple {viV2---v,n) where 
^ ^ Vi < di, i = 1,2, ■ ■ ■ ,171. In other words, we set up (G, a) using the standard basis in 
H. Since (G, a) is the graph of a pure state, it must be a clique on some subset of V{G, a), 
all vertices not belonging to this subset being isolated [2] . We divide the m parts of the 
system in two nonempty disjoint subsets (partitions) whose union makes up the whole 
system. We call them s and t, where t is the complement of s in the set of all parts of the 
system. That is, s and t are the nonempty subsets of {1, 2, ■ ■ ■ , m}, sUt = {1,2, ■ ■ ■ , m} 
and s n t = 0. This corresponds to H = H^'^ ® n^^\ H^'^ = H"'^ (g) W^'^ ® U'^Ki o • • • o 
n< , T^W = Ti'^h n'^j2 ® Ti'^js • • • Ti'^h , {ii, i2, • • • , is} = s, {ji,j2, • ■ ■ , it} = t. As a 
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result of this division, we can divide the m tuple v = {viV2 ■ ■ ■ Vm) into the corresponding 
partitions (strings) which we call Vs,Vt. Thus Vs — Vi^Vi^ ■ ■ - Vi^ and Vt — vj^Vj^ • • • Vj^. We 
can label each vertex equivalently by {vs,Vt). We call Vs The s part and Vt the t part 
of the vertex label. For example, consider a four partite system whose parts are labeled 
1, 2, 3, 4. Let s = {1, 4}, t = {2, 3}. Then a vertex label (1122) can be written as (12, 12). 
A vertex label {V1V2V3V4) becomes (f 1^4, ^2^3)- 

Partial transpose with partition s. This is a graph operator denoted Tg operating 
on E{G, a) which we define separately for the graphs with real and complex weights. 

Definition 2.1 : Let (G, a) be a graph with real weights. The operator Tg is defined 
as follows. 

T, :(G',a)i — >{G^%a'), {G,a) 3 {{vs,Vt),{ws,Wt)} 1 — > {{ws,Vt), {vs,Wt)} G (G'^^a'), 
with 

a'{{{ws, Vt), {vs, wt)}) = a{{{vs, Vt), {w^, Wt)}), 

that is, a'{Tse) = a(e). Note that, in general, E{G, a) is not closed under T^. The operator 
Tt for the partition t giving the complement of s in the m-partite system is defined in the 
same way. Note that Tg = Tt. 

Definition 2.2 : Let {G, a) be a graph with complex weights. The operator is 
defined as follows. 

T, :(G',a)i — > {G^%a'), (G, a) 3 {{vs,Vt),{ws,Wt)} 1 — > {{ws,Vt), {vs,Wt)} G (G^^a'), 
with 

\a'{{{ws,vt),{vs,wt)})\ = \a{{{vs,vt),{ws,wt)})\, 

that is, |a'(T5e)| = |a(e)|. Note that, in general, E{G,a) is not closed under Tg. The 
operator Tf for the partition t giving the complement of s in the m-partite system is 
defined in the same way. Note that Tg = Tt. 

Definition 2.3 : Partial transpose of (G, a) with real or complex weights with respect 
to partition s, denoted {G^% a') is the graph obtained by acting Tg on E{G, a). Note that 
V{G, a) = V{G^%a'), but in general E{G, a) ^ E{G^%a'). 

An edge joining the vertices, whose labels have either the same s part or the same t 
part or both are fixed points of Tg. We have 

Tg{{vg,vt),{vg,wt)} = {{vg,vt),{vg,wt)}, 

with the condition on weights automatically satisfied. Similarly 

Tg{e) = Tg{{vg,vt),{wg,vt)} = {{wg,vt),{vg,vt)} = {{vg,Vt), {wg,Vt)}, 

with the condition on weights automatically satisfied. Note that in the case of complex 
weighted graphs the action of Tg on {{vg,Vt), {wg,Vt)} changes its orientation and hence 
a'{Tg(e)) = a*(e). But by the definition of Tg this edge is still invariant under Tg. 
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Further, the definition of the operator Ts leaves the phase of a'(Ts(e)), e G E{G,a), as 
a free parameter. We shall use this freedom later in fixing the phases of the graphs 
corresponding to the factors in the tensor product decomposition of a density matrix. 

If both s and t parts of two vertices of e G E{G, a) are the same then both vertices 
are identical and wc have a loop, which is obviously preserved under T,. Thus divides 
E{G, a) into two partitions, one containing all fixed points of Tg, that is, edges with same 
s or t part and all loops, which we call T set and the other containing the remaining edges 
which we call C set. In other words {G, a) is the disjoint edge union of the two spanning 
subgraphs corresponding to the set and the C set. 

Note that is the identity operator {T^ = 1). Thus Tg is its own inverse and is one 
to one and onto. 

Lemma 2.1 : Let {G, a) be a graph of a pure state in the Hilhert space H — Ti.^^ (8) 
7^'^^ (g) • • • (g)?^'^™'. Let {G'^% a') be the partial transpose of {G, a) with respect to a partition 
s as defined above. Then E{G, a) is closed under Tg, E{G, a) — E{G'^% a'), if and only if 
A(G',a) =^{G^%a'). 

We emphasize that the closure of E{G^ a) under T, means, for every e G E{G, a), Ts{e) = 
e' G E{G,a) and a(e') = a(e) or |a(e')| = |a(e)| as appropriate. 

Proof : Only if part : We are given that E{G, a) is closed under Tg. We divide E{G, a) 
into two partitions C and as above. Note that if E{G, a) is closed under Tg then the 
sets C and are separately closed under Tg. Consider now the set of edges incident on 
a vertex in V{G,a). The edges in this set which belong to J-' set are not shifted by Tg. 
Since E(G,a) is closed under Tg and — 1, every incident edge belonging to C is the 
image of an edge in C with same weight ( or same absolute value for the weight) under 
the action of Tg on C. Thus the degree of each vertex is preserved under the action of Tg 
on E{G, a), for both, real and complex weighted {G, a) [2], so that A{G, a) = A{G'^% a'). 

If part : We are given A{G,a) = A{G^%a'). The edges in E{G,a) belonging to T 
remain in E{G,a) under the action of Tg. Now suppose that the set C is not closed under 
Tg. Since Tg is its own inverse, the edge e for which T's(e) ^ E(G, a) cannot be the image 
of any other edge in E{G,a) under Tg. Therefore, the degree of the end vertices of e is 
changed under the action of Tg. This contradicts the assumption A{G,a) — A(G^-', a').M 

Lemma 2.2 : Let {G, a) be the graph of a pure state in the Hilbert space H — Tif^^ ® 
a m-partite quantum system. Then {G,a) — {Gg,b) □ (Gt,c) where 

(Gg, b) is the graph of a pure state in the Hilbert space made up of s factors 'H^^'' — 
®l-L'^'-2 (g)7Y'''3 (g) • ■ ■ (g)7i'^'s and {Gt, c) is the graph of a pure state in the Hilbert space 
made up of t = m — s factors Tif-^^ = Ti'^-'i ® Vf'^'^ ® ® ■ ■ ■ ® 7^°'^* , if and only if the 
edge set E{G, a) is closed under Tg, s — {ii, i2, is, • " " > 'is}- 

Proof : Only if part : Case I : Graphs with real weights. We are given 

{G,a)^{Gg,b)B{Gt,c) 
= C{Gg, b) ^Cr]{Gt, c)+C{Gg, b) ®N{Gt, c) b) ®£(Gt, c) b)m{Gu c) (1) 
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Using the definition of the operators C,ri,J\f,fl [2] and that of the tensor product of 
graphs, we can make the following observations. The second term is a spanning subgraph 

of {G, a) each of whose edges has a common t part and hence is fixed point of Tg. The 
third term is a spanning subgraph of {G, a) each of whose edges has a common s part, so 
that each edge is a fixed point of T,. The fourth term is a spanning subgraph of (G, a) 
which contains only loops all of which are fixed points of Tg. Again, from the definition 
of the tensor product we see that in the first term, any edge {vs,Ws} in jC{Gs,b) with 
weight b{{vs,Ws}) and any {vt,Wt} in Crj{Gt,c) with weight c{{vt,Wt\) gives us, under 
the tensor product, two edges {(^s, ft), (w^, Wt)} and {(fs, Wf), (tt;^, ff)} with the same 
weight h{{vs,Ws})c{{vt,Wt}), which are the images of each other under Tg. This proves 
that E{G, a) is closed under Tg. 

Case II : Graphs with complex weights . We are given 

{G,a)^{Gs,b)B{Gt,c) 

= C{Gs, b) ® C{Gt, c) + C{Gs, b) ® M{Gt, c) + M{Gs, b) ® C{Gt, c) 

+MGs, b) ® n{Gt, c) U 2J\fjr{Gs, b) ® AfCrjiGt, c)} (2) 

Note that the first three terms are similar to those in Eq. (1) and the arguments corre- 
sponding to these terms in the paragraph following Eq.(l) apply, except that we require 
|a(Tse)| = |Q(e)|. Again, fourth term correspond to graph with loops (and no edges) which 
are fixed points of Tg. This proves that E{G, a) is closed under Tg. 

If part : We begin by noting that the graph {G, a) has the structure of a clique 
and isolated vertices, with a)| = d = did2---dm where di is the dimension of 

H'^S 1 < i < m, because it is the graph of a pure state in 7i = H'^^ ®'H'^'^ ® ■ ■ ■®H'^'^. Let 
{Kn, a) denote the clique in [G, a) and 14 (G, a) be the set of vertices on the chque. Let 
|Vfc(G,a)| = n. We are given that E[G,a) is closed under Tg. Note that all loops arc on 
the clique and no loops are on the isolated vertices. Consider a vertex {vg,Vt) on {Kn,a). 
Let q denote the number of vertices in {Kn,a) having the same s part as {vs,Vt) and p 
denote the number of vertices on {Kn, a) with the same t part as {vg, Vt). We note that p 
and q are the same for all vertices on [Kn, a), otherwise the set C is not closed under Tg. 
We draw {K^ a) as a lattice of p rows and q columns, such that all vertices in one row 
have common s part and all vertices in one column have common t part. Since {K^. a) 
is a complete graph, from figure 1, we see that any vertex {vg,Vt) has {p — ^) + {q — 1) 
neighbors giving edges in the T set and {p —l){q — 1) neighbors giving edges in the C set. 
Since (/^n, a) is complete (v^, Vt) has n — 1 neighbors giving n = pq. 

Now consider C set on [K^, a). Prom the definition of the tensor product of weighted 
graphs [2], we can factorizc each pair {{vg,Vt){tUs^Wt)}\ {(ws,ft), (■^s,'"^*)} in the set C as 
the tensor product of two edges {vg, w^} and {vt^Wt} with weights b' and d satisfying 

a{{{vg, Vt), {wg, wt)}) = b'{{vg, Wg}) ■ c'{{vt, Wt}) = a{{{wg, Vt), {vg, Wt)}) 
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figure 1: Every row contains vertices with same s part, and every column contains 
vertices with same t part. In a vertex label, first number stands for the s part and the 

second for the t part. For example, the edges between the vertex 34 and all the vertices 
in the 3rd row and 4th column are in the set JF while the edges between vertex 34 and 
all the vertices in the four blocks obtained by deleting 3rd row and Ath column are in 

the set C. 



or 

\a{{{vs,vt),{ws,wt)})\ = \b'{{vs,Ws})\ ■ \c{{vt,wt})\ = \a{{{ws,Vt), {vs,Wt)})\. 

Writing each pair {e, T^e} in C set in this way and taking the disjoint edge union [2] in 
all these tensor products we get {Gc, a) — {G\, h') ® {G'^, c') where (Gc, a) is the spanning 
subgraph on {Kn,a) corresponding to set C. {G'i,h') and {G^ic') are graphs on p and 
q vertices respectively. Again, from the definition of the tensor product of weighted 
graphs, we know that isolated vertices in the factors produce isolated vertices in the 
product. Therefore we can add [dg — p) isolated vertices to [G'^, h') and [dt — q) isolated 
vertices to (G"2, c'), where dg — di^di^ ' ' ' dis ^^id dt — dj^dj,^ ■ ■ -dj^. We call this graphs 
C{G'g,b') and CiG^^d) (the operator C defined in [2] removes loops from a graph). The 
tensor product C{G'g,b') (E) C{G[,c') gives the spanning subgraph of {G,a) corresponding 
to set C. Now consider a row in figure 1 containing vertices with common s part say 
Vg- This row generates q{q — l)/2 edges of the form {{vs,Vt), {vs,Wt)} all in the set J^. 
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By the definition of tfie cartesian product of the weighted graphs [2] , each of these edges 
is the cartesian product of the vertex Vg in say {Gi,b) with the edge {vt,Wt}, in say 
(G2, c) where a{{{vs, Vf), {vs, Wt)}) = c)„^c({vt, Wt}). Thus the graph {G2, c) is a graph of q 
vertices obtained by projecting each of the q{q — l)/2 edges {{vg, Vt), {vg, Wt)} to {vt, Wt}, 
with corresponding weight assignments, and thus is a complete graph on q vertices. As we 
vary Vs through its p possible values, one row corresponding to each value, the definition of 
the cartesian product of weighted graphs generates the same ((^2, c) possibly with different 
weights on edges. In exactly the same way, q columns in figure 1 generate the complete 
graph on p vertices (Gi, b) by employing the cartesian product of weighted graphs. Noting 
that cartesian product of an isolated vertex in (Gi, b) with an edge in (6*2 , c) or vice versa 
gives an isolated vertex in the product graph, we have shown that the spanning subgraph 
of (G, a) corresponding to the set edges gets generated by &)n(Gt, c) containing 
ds and dt vertices respectively. Finally note that there are no loops in {Gs,b)n{Gt,c) 
because loops contribute to the cartesian product of weighted graphs only via the degrees 
and 9^,4. In the above analysis, question can be raised regarding the weight functions 
of the factors. Thus more than one weight functions can generate the spanning subgraph 
corresponding to the C set while the T set edges on different rows in figure 1 may be 
generated by different weight functions on the factors in the cartesian product. These 
points will be addressed later in this proof. We denote by Gi and G'l to be the graphs 
underlying {Gi,b) and {G[.b') and similarly define G2, G2, G, G[ ® G2 and Gi\I\G2. 
We first note that \V{Gi)\ = p = \V{G[)\ and \V{G2)\ = q = \V{G'2)\. In fact, we 
can identify the set of vertices for Gi and G[ and for G2 and G'2 respectively, that is 
V{Gi) = V{G[) and ^(Ga) = ^(G's). We now show that the identity map on V{Gi) is 
an automorphism taking Gi to G'l. In other words, Gi and G'l are identical. Consider 
a vertex {Vs,Vt) in G and let NQ{vs,Vt) denote its neighborhood in G. We denote by 
NGinG2{vs,Vt) and NG'^^G'^{vs,Vt) the neighborhoods of {vs,Vt) in GinG2 and G[ G2 
respectively. In other words NG[(s,G'^{vs,Vt) contains neighbors of {vs,Vt) with edges in 
set C and NG^nG2i''^s,Vt) contains the neighbors of {vs,Vt) with edges in set J^. Clearly 
A^G'i®G^(^s,^'t) and NG,nG2i'"s,Vt) partition Naiv^Vt). Let VxiGi) and Vk{G2) be the set 
of vertices on the clique in Gi and G2 respectively. Consider [14] 

NG[^G'^{Vs,Vt) = Va-(G) - {iVs,Vt)} - NG,nG2{Vs,Vt) 
= Vk{G) - {Vs} X {vt} - {{Vs} X NG^iVt) U Ng,{Vs) X {vt}} 

= Vk{G,) X Vk{G2) - K) X {vt} - {{vs} X NgM U NgM x {vt}} 

= {Ng, (Vs) U {Vs}} X {Ng2 (vt) U {vt}} - {Vs} X {vt} - {{Vs} X Ng2 (Vt) U Ng, (Vs) X {vt}} 
= NgAVs) X NgM U NgM X {vt} U {Vs} X NgM U K} X {vt} - K) X {vt} 
-{{Vs} X NG^ivt) U NgM X {vt}} 
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Therefore [14], Ng'^{vs) x No'^ivt) = x 

or NG'^ivs) = NG.ivs) and Na^ivt) = NgM)- 

Therefore, Gi and G[ and G2 and G2 are identical. In particular Gi and G2 are 
cliques. Gg — G'^ and = G[ are chques plus isolated vertices. Thus we get 

(G, a) = 6') c') + (G„ h)n{Gu c) + (G", a) (3) 

where (G", a) is the graph obtained from (G, a) by removing all edges and keeping loops. 

The only remaining gap is to show that a consistent assignment of weights to the 
factors Gs and Gt is possible so as to express (G, a) as modified tensor product. To 
get the required weight assignments we use the requirement that both the factors in the 
modified tensor product must correspond to pure states. Indeed, we know that both Gg 
and Gt have the form of clique plus isolated vertices as required for them to represent 
pure states. We know that the graph (G, a) corresponds to pure state. Therefore, its 
weight function a must satisfy [2], assuming {G,a) to be a real weighted graph, 

(.Vs,vt)eV{G,a) eeE(G,a)) 

Splitting E{G, a) into C and T sets and using the definitions of the tensor and cartesian 
products of weighted graphs, we get, (note that a paired label is for a vertex in V{G, a) 
and single labels with suffix s and t are vertices in V{Gs.i b) and V^(Gt, c) respectively), 

E ^(^..'^t) + 2 E E {a^{{{vs,vt),{ws,wt)}) + a^{{{vs,wt),{ws,vt)})) 

{Vs,Vt) (.Vs,Vt) {Ws,Wt) 

Ws^Vs,Wt^Vt 

+2 E ^ a\{{vs,vt),{vs,wt)}) + 2Y^ ^ a\{{vs,vt),{ws,vt)}) = i>lG,ay 

{Vs,Vt) {Vs,'Wt) (.Vs,Vt) {Ws,Vt) 

Since E{G, a) is closed under Tg we can write 

E ^%s,vt) + 4 E E "'(^ ^*)' ^*)» 

(Vs,Vt) (Vs,Vt) (Ws,Wt) 

+2 E E E E = (5) 

(Vs,Vt) (Vs,Wt) (Vs,Vt) (Ws,Vt) 

Using splitting of the weight functions in C set and T set we get 
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+2E<E E c^(K«^*}+2E<E E ^'(K«^J) = ^ka)- (6) 



lis Vt Vt Va Ws 



Since the graphs (G^, b), {Gt, c), (G^) b'), and (G^, c') correspond to pure states, b', c', b, c 
must satisfy, 

E< + 2E E b\{vs,w.})-X>laa,) (7) 
E<+2E E b'\{vs,w,})=^la^,,^ (8) 

E<+2E E ^'(K^*})=^kc) (9) 



Vt Vt 

Wt^Vt 



E<+2E E = ^(G.co- (10) 



Wt^Vt 



We see that Eqs. (7), (8), (9), and (10) are consistent with Eq. (6) provided 
(0 ^{vs,vt) = '^vs^vt for all {vs,vt) e V{G,a) and consequently D(G,a) = 3(G.,6)f(G(,c) 
(ii) = c^{{vt,Wt}) = c''^{{vt,Wt}). 

We first fix a vertex {vs,Vt) G V{G,a) and obtain its degree <){vs,vt)- Summing the 
edges with the same s part we get 

a{{{vs,vt),{vs,wt)}) = 'Ovs E c{{vt,wt}). (11) 

{Vs,Wt) ""J 

Wt^Vt 

Adding edges with the same t part we have 

E ^{{{'"s,vt),{ws,vt)]) = -0^^ ^ 6({ws,wJ). (12) 

Adding over edges in the C set we get 

a({(^;„'f;t),(ws,'u;t)}) = b'{{vs,Ws}) ^ c'({T;t, Wt}). (13) 



Adding these three terms and the weight of the loop on {vs,Vt) we get (){vs,vt)- 

The 

requirement that ()(vs.vt) = ^v^^vt is satisfied provided 

[iii) b'{{vs,Ws}) — b{{vs,Ws}) and c'{{vt,Wt}) — —c{{vt,Wt}) which leads to 

{■iv) ^{vs,vt) - a{{{vs,vt), {Vs,Vt)}) = Hv.Vy^ - b{{Vs,Vs})c{{Vt,Vt}). 

This is satisfied provided a{{{vs,Vt), {vs,Vt)}) = b{{vs,Vs}) ■ c{{vt,Vt}). This require- 
ment is consistent with {G", a) — fl{Gs, b) (8) fl{Gt, c). 
We can write, therefore, 

(G, a) = LiGs, b) LiGt, -c) + (G,, b)D{Gt, c) + Q(G'„ b) Q(G't, c) 

^{Gs,b)B{Gt,c) 

Now, let us deal with the case where {G, a) is a graph with complex weights [2] . In 
this case Eq. (4) is replaced by 

(.Vs,vt)eV{G,a) eeE(G,a)) 

and Eqs. (5, 6, 7, 8, 9, 10) become 

E E E H{{Vs,Vt),{Ws,Wt)})\^ 

Ws^Vs,Wt^Vt 

+2 E E H{(Vs,Vt),(v,,Wt)})\' + 2j2 E l«({K,^t),K,^t)})r = f(G,a) (5&) 

{Va,Vt) (vs,'Wt) {Vs,Vt) (Ws-Vt) 

Wt^Vt Wa^Va 



E 4.M+4E E i^'(K,«^J)i'E E 

{Va,Vt) Va Vt 

+2E<E E ic(K^*})r+2E<E E iKK^J)r=^(G,a) (6^) 

Vs Vt V>t Wa 

Wt^Vt Ws^Vs 



E<+2E E iKK^J)i'=^k^» i^b) 

Vs Vs 

WsJ^Vs 
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Vs Vs 



Vt Vt 

Wt^Vt 



E<+2E E i^'(K'^*})r=fkc')- (10^) 



Vt Vt 

wti^vt 



We see that Eqs. (7b), (8b), (9b), and (10b) are consistent with Eq. (6b) provided 

,Vt) ^Vs^Vt 

for all {vs,vt) e V{G,a) and consequently c)(G,a) = '^{Gs,b)^{Gt,c) 
{vi)\h{{vs,WsW = \b'i{vs,w,})\'; |c(K«;,})p = \c'{{vt,Wt})\\ 
Eqs. (11, 12, 13) become 

E klil^s,-"*), (^^s,wt)})| = c)„, ^ |c({vt,wt})| (116) 

H{{vs,vt),{ws,vt)})\ = 0., J] |6(K,«'J)| (126) 



J] |a({(^;„^;,),(^/;„ «;,)})!= ^ \h' {{v,,Ws})\ E |c'(K,^t})l- (136) 

(l's,1ft) 
Wt^Vt 



Wt 



Adding these three terms and the weight of the loop on {vs,Vt) we get ^{vs,vt)- The 
requirement that '^{ys,vt) — '^vs^vt is satisfied provided 

{vii) \h'{{vs,Ws})\ = \b{{vs,Ws})\ &nd \c'{{vt,wt})\ = |c({^;^, Wt})| 
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{via) a{{{vs,vt),{vs,vt)}) = b{{vs,Vs})c{{vt,Vt}) - 2 ^ \b{{vs,Ws})\ ^ \c{{vt,Wt})\. 

Ws Wt 
Ws^Vs Wt^Vt 

The requirement {viii) is consistent with 

{G\a) = n{Gs,h)®n{Guc)VA2MC{Gs,h)®MCr]{Gt,c). (14) 

We now choose phases of the weight functions b and c. Consider the edges e = 
{{vs,vt),{ws,wt)} and T^e = {{ws,Vt),{vs,Wt)] in E{G,a). We know that |a(e)| = 
|a(T5e)| . Let 6*^^ and e*^^ be the phases of Oi = a(e) and a2 = a,{Tse) respectively. If we 
require that these two edges in E{G, a) be produced by the tensor product of the edge 
{vg, Ws} in {Gg, b) with the edge {vt, Wt} in {Gt, c) (see figure 2) then the phases of weights 
b and c on the corresponding edges must be 0i = (^i + 62) /2 and 02 = (^1 — 02)/2 
respectively. 




figure 2 

This completely fixes the weight functions b and c on {Gg, b) and (Gj, c) respectively. 
We now have, for every edge e in E{G,a) the corresponding edges 61,62 in {Gs,b) and 
(G(,c) respectively such that 6 = 61 (g) 62 and a(6) = &(6i)c(62). Thus we have, finally 
from Eq. (3), {vii) and (14) 

{G,a) = C{Gs,b)®C{Gt,c)^C{Gs,b)®N{Gt,c)+N{Gs,b)®C{Gt,c) 

H^{Gs, b) ® n{Gt, c) U 2AfC(Gg, b) (g) MCr,{Gt, c)} (2) 
= (a,6)H(G„c). 

■ 

Lemma 2.3 : Let o", and at he density matrices for pure states. Then a = (jg® at 
if and only if {G, a) = {Gg, b) □ {Gt, c) where {G, a), {Gg, b) and {Gt, c) are the graphs for 
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cr, Gs and at respectively. 

Proof : This Lemma is identical with theorems 4.5, 6.10 in [2]. 



Theorem 2.1 : Let (G, a) be the graph of a m-partite pure state a in the Hilhert space 
'H'^^®H'^^®---®n'^^. Letn^'^ =7i'^n®7i"'«2®7i"'«3®...®7i*., {ii,i2,--- ,ij = s cor- 
respond to an s set, H^*^ — T^'^-'i ®1-L'^^'^ ^l-C^^^ ® ■ ■ ■®'H'^^t , {ji, j2, • • • ,3t\ — t, correspond 
to the t set which is the complement of s set in {1, 2, ■ ■ ■ , m}. Then a = a^® at, where 
as and at are pure states in H^^'' and liS^'> with graphs {Gg, h) and {Gt, c) respectively, if 
and only ifA{G,a) = A{G'^%a'). 

Proof : Using Lemmas 2.1, 2.2 and 2.3 we have 

a = as ^ at <^=^ (G, a) = (Gg, b) □ {Gt, c) <^==^ EiG, a) is closed under Tg <^=^ 
A{G, a) = A{G'^%a'). A state a is entangled if A{G, a) ^ A(G^'', a') in every partition s 
and t of {1, 2, • • • , m}. ■ 



III. ALGORITHM 



While proving the if part of Lemma 2.2, we have shown that the number of vertices 
in the cliques of the factors Gi and G2 defined there {p and q respectively) are the 
factors of the number of vertices on the clique in (G, a), that is, n = pq. This means 
that the m-partite pure state j'^) corresponding to (G, a) has two factors and \ip2), 
corresponding to Gi and G2 respectively, such that lipi) lives in a p-dimensional subspace 
of 7i*i C?>?i*3 ®- ■ and \^2) lives in a g-dimensional subspace of T-C^n ^l-f'-h 

T-C^^i (8) ■ ■ ■ ® 7i"'j« . If the weighted versions of Gi and G2, namely (G^, b) and {Gt, c), can 
be further factorized, the dimensions of the corresponding subspaces will be the factors of 
p and q respectively. This procedure can be iterated at most until the dimensions of the 
subspaces for the factors of are the prime factors of n. Therefore, the dimension of the 
subspaces containing the factors of {ip) arc the prime factors of n or the products of such 
factors. This fact can be used to get a polynomial algorithm to find the full separability 
of 1'^) in the following way. By full separability we mean expressing as a product 
state whose further factorization is impossible. Denote hy pi > p2 > ■ ■ ■ > Pk the prime 
factors of n. Let di > d2 > ds > ■ ■ ■ > dm he the dimensions of the Hilbert spaces of 
m parts arranged in a non-increasing order. Let Si > be the least integer satisfying 
did2 - ■ ■ ds^ > Pi- We implement our algorithm (Theorem 2.1) on partitions {s,t) with 
Si < s < m — 1. The total number of times the algorithm has to run, in the worst case, is 




which is a polynomial of degree Si in m. Thus we have a polynomial algorithm to 
check separability of a m-partite system. Suppose we get the separability as H^^^ (8> 7Y^*^ . 
Then the factor in Ti^^^ cannot be further factorized as it corresponds to the largest prime 
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factor of n and 7i*^*' contains factors corresponding to p2 > Pa • • • > Pfc- We repeat the 
above algorithm on TC^^^ with p2 as the largest prime factor. Its worst case complexity is 
given by a polynomial of degree (m — s)*^ where S2 is defined like Si above. Wc carry out 
these iterations until full separability is obtained. Thus if we do not get any factorization 
in the first iteration, corresponding to the largest prime factor pi, then the state is fully 
entangled, like GHZ or W state. Unless the factorization carries up to m factors, the 
factors of the state contain one or more entangled states involving less than m parts. The 
total algorithm is polynomial in m. Note that, if n is prime, then all that is necessary 
is to look for some Vi common to the m tuples for all vertices on the clique. If say Vi is 
common, then 

\^) = \<P)®\Vi) 

with 10) eW^' W^'-^ ® H'^'+^ (8) • • • ® W^"" and \vi) e H'^\ Otherwise the given 
state \ip) is entangled. 

IV. A COUNTER EXAMPLE 

We note that Theorem 2.1 may not apply to mixed states as the following exam- 
ple shows. Consider the bipartite separable state a — l/2\y, —)\y, —){y, —\{y, —\ + 
l/2\x,+)\x,+){x,+\{x,+\, where \y , -) = j^{\0) - and |x, +) = ^(|0) + |1)). The 
corresponding density matrix in standard basis is 

2 1 + i 1 + i " 
1-i 2 2 1+i 

l-i 2 2 1+i 
l-i l-i 2 

and the corresponding graph is shown in figure 3. We see that the C set contains only 

one edge for all possible partitions and hence cannot be closed under any Tg. We show, in 
section (V), that the degree conjecture applies to states with real weighted graphs without 
loops. Therefore, the above example shows that the degree conjecture does not apply to 
all mixed states. 

V. PROOF OF DEGREE CRITERION FOR A CLASS OF MULTIPAR- 
TITE MIXED STATES 

In this section we prove the degree criterion for the class of states whose graphs have 
no loops and have real weights. 

Theorem 5.1 : Let (t{G, a) be a density matrix acting on H.'^^ ®1-L'^'^ ® ■ ■ with 
a real weighted graph (G, a), ond = ^1^2 ■ ■ - dm vertices, having no loops. If a) is sep- 
arable in s, t cut where s, t is a partition of {1,2,..., m}, so that cr(G, a) = Picr^^^ ®<^f\ 
where af'' and af^ are density matrices acting on Ti^^^ and 7i^*^ with graphs {G^i \ b) and 
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figure 3 

{Gf\c) respectively, then A(G', a) = A(G'^^ a'). 

Proof : Let Q{G, a) be the Laplacian of a graph {G, a) with real weights without 
loops, on d vertices. For a graph without loops Q{G,a) = L{G,a) [2]. Let D be any 
d X d real diagonal matrix in the standard orthonormal basis i — 1,2, . . . ,d, such 

that D ^ and Tr(D) = 0, where Tr is the trace of D. This means that there is at 
least one negative entry in the diagonal of D. Denote this element by Da = 6j < 0. Let 
l^o) = Ej ki) and 10) = T.jXj\vj) where, 

_ r iij^z 

\ k ii j = i, k real. 
Let \x) = iV'o) + 10) = E-=i(l Then 

{x\L{G,a) + D\x) = {^o\LiG,aMo) + {HL{G,am + {c^\LiG,a)\^o) + 

{4>\L{G,aM + (^ol^l^o) + {^oW) + im^o) + (01^10) (15) 

Since L{G, a) is positive semidefinite, we must have, for the quadratic form associated 
with L{G, a) of a graph without loops [15,2], 

x^L{G,a)x — ^jii^j ~ ^iT ^ 0) (16) 

{j,l}€E{G,a) 

for every x G M'^. Since \iIjq) is (unnormalized) vector having all components equal to 
unity, from Eq.(16) it follows that (^o|^(G', a)|V'o) = 0. Also (V^ol^l^o) = Tr{D) = 0. 
We have 

{cl>\L{G,aM^e{L{G,a))u^k\, 

where tJj denotes the degree of ith vertex. For a real weighted graph without loops the 
sum of the elements in any row of its Laplacian is zero. This leads to 

{iJo\L{G,a)\<P) = (0|L(G,a)'^|^o) = (0|L(G', a)|^o) = 0. 
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Finally, the remaining terms in Eq. (15) are given by 

{ct>\D\ct>) = he 

{MD\^) = hk = {ct>\D\ij^). 

Thus 

{x\L{G,a) + D\x) = k\h + l),) + 2khi 
So we can then always choose a positive /c, such that 

(x|L(G',a) + L>|x)<0. 

Then it follows that 

L(G,a)+D^O. (17) 

This expression is identical with that obtained in [1,2]. For any graph (G, a) on 
d — dgdt vertices 

consider the degree condition A(G, a) = A{G'^%a'). Since M'^'{G, a) = M{G'^%a') where 
M is the adjacency matrix, we get, 

(L(G, a)f^ = {A{G, a) - A{G^%a')) + L{G^%a'). (18) 

Let 

D^A{G,a)-A{G^%a'). 
Then D is an d x d real diagonal matrix with respect to the orthonormal basis 

\vi) = \ui) (g) . . . , \vd) = \ud,) (8) \wdt). 

As the degree sum of a graph is invariant under partial transpose, 

Tr{D) = Tr{A{G, a)) - Tr{A{G^%a')) = 0. 

We have two possible cases : D ^ or D — 0. If D 7^ 0, that is, the degree condition 
is not satisfied {i.e.A{G,a) ^ A{G'^%a')), we can write, via Eq.(17), L{G'^%a') + D ^ 
because, {G^%a') is a real weighted graph without loops and D 7^ with Tr{D) = 0. 
Using Eq.(18), this means {L{G,a)f^ %_ 0. As a(G, a) = ^L(G',a) [2], (T{G,a) is 

entangled [5]. 

■ 

In order to test the separability of the m-partite state o"(G,a), we have to apply the 
degree criterion to all the bipartite cuts (s, t) of m-partite system. This procedure may 
not detect the full separability of m-partite state [16]. However, all the known separability 
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criteria for multipartite states test separability only in bipartite cuts, and hence are not 
enough to guarantee full separability. 

In this paper we have settled the degree conjecture for the separability of multipartite 
quantum states. Recently, Hildebrand, Mancini and Severini have proved that the degree 
criterion is equivalent to the PPT criterion for separability [17]. However, the importance 
of degree conjecture ensues from the opportunity it offers to test the strengths and limi- 
tations of a nascent approach to the separability problem. We see that this approach has 
contributed to test the separability of a class of multipartite mixed states (Theorem 5.1) 
as well as to the efficient factorization of multipartite pure states. 
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